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We provide numerical evidence that the low-lying part of the entanglement spectrum of a real- space block 
(i.e. a single interval) of a one-dimensional quantum many body system at a conformal critical point corresponds 
to the energy spectrum of a boundary conformal field theory (CFT). This correspondence allows to uncover a 
subset of the operator content of a conformal field theory by inspection of the entanglement spectrum of a single 
wave function, thus providing important information on a CFT beyond its central charge. As a practical applica- 
tion we show that for many systems described by a compactified boson CFT, one can infer the compactification 
radius (governing e.g. the power law decay of correlation functions) of the theory in a simple way based on the 
entanglement spectrum. 
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Introduction. — The study of entanglement properties of 
quantum many body systems has established itself as a pow- 
erful approach to uncover many new aspects of strongly cor- 
related quantum systems [1,2]. An initial focus was on the 
understanding of entanglement entropies, but in recent years 
the concept of the entanglement spectrum (ES) - introduced 
by Li and Haldane in Ref. [3] - attracted also significant in- 
terest. One of the main questions asked in this context is 
whether the spectrum of a reduced density matrix contains 
more physical information than an entanglement entropy. For 
fractional quantum Hall wave functions for example, the en- 
tanglement spectrum contains additional information on the 
theory describing the chiral edge modes, which is not acces- 
sible in entanglement entropies [3]. In the meantime the ES 
has become a powerful tool to analyze topological states of 
matter and symmetry protected topological phases (such as 
topological insulators and Haldane insulators) as well as non- 
topological states, such as e.g. continuous symmetry breaking 
states [4]. 

One of the first areas where entanglement entropies were 
calculated, are systems described by conformal field theory, 
especially in 1+1 dimensions, where the conformal symme- 
try is particulary rich [5-7]. The properties of entanglement 
entropies in 1+1 dimensional CFTs are quite well understood 
by now [8-16], and perhaps best highlighted by the celebrated 
formula for the entanglement entropy between an interval A 
and its complement in a finite system [here given for the case 
of periodic boundary conditions (PBC)]: 
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where pA = T^^L\A\i^) (^|. {L) is the length of the interval 
A (the entire system), and s[ is a non-universal constant. An 
important feature of this formula is that it only depends on the 
central charge c, and no other property of the CFT enters. 

In the context of 1+lD CFT it is still an open question 
whether the entanglement spectrum {^}, i.e. minus the log- 
arithm of the eigenvalues of {A} = {exp(— ^)} con- 
tains additional useful physical information beyond the cen- 
tral charge. Calabrese and Lefevre [17] calculated the density 
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Figure 1. (Color online). Setup of the bipartitions and boundary 
conditions of the one dimensional lattice models investigated here. 
The orange shaded region denotes the size and position of block A 
in the case of (a) open boundary conditions (OBC) and (b) periodic 
boundary conditions (PBC). 



of states of the entanglement spectrum in CFT, which turned 
out to be universal and to depend only on the central charge. 
So the answer seems to be negative at first sight, but their re- 
sults applies to the inner, continuous part of the entanglement 
spectrum, and not directly to the discrete part at the lower end 
of the spectrum. 

It is the aim of this paper to provide compelling numerical 
evidence that the lower part of real-space ES contains inter- 
esting CFT structure by showing several examples where the 
ES can be understood as the energy spectrum of a boundary 
CFT on a strip of effective width ^ logL. Since the energy 
spectra of boundary CFTs are organized in terms of conformal 
families [18-20], one obtains access to a subset of the oper- 
ator content of the underlying CFT by inspection of the ES. 
In the case of c = 1 theories corresponding to a compactified 
boson, a recurrent task is to infer the compactification radius 
of the theory for a given microscopic model. Here we show 
that the boundary conditions of the boundary CFT and there- 
fore the operator content of many interesting systems are such 
that the compactification radius (governing e.g. the exponent 
of correlation functions) can actually be deduced by a simple 
procedure from the ES. 

We start by describing the simulation and entanglement 
spectrum setup and then discuss the observed operator con- 
tent of the entanglement spectrum in various Luttinger liquids 
such as the Bose-Hubbard chain, the S = 1/2 XXZ chain 
and the 5* = 1 Heisenberg chain in a magnetic field, as well 
as the the transverse field Ising model and the quantum three- 



2 




2 



O U/J=2, L=256 Obc 



(b) (c) 



" JJLP 6 



-3-2-10123 

5N, 




sJLP 

< 

3JLP 

3JLP 

^^-Ot<^^-0-^-O- 1 

O U/J=2, L=256 pbc 



-3-2-10123 

8N, 



OS0S0SO 
0808080 

0§0§0 

0g0 , 



-3-2-10 1 2 3 
(5A^A or (^S'^ 



Figure 2. (Color online) (a) / (c) Entanglement spectrum ^ for a L = 256 Bose-Hubbard chain at [// J = 2 with OBC / PBC and a block size of 
La = 128. (b) / (d): Structure of the ES in each 6Na sector for OBC / PBC, obtained after subtracting the lowest ^ value in each 6Na sector, 
and setting the difference between the first and second ^ value in the 6Na = sector to one. The levels with an assigned value are indicated by 
the dashed red line, while the relative position of all other levels highlights the emergent CFT structure of the ES. The approximate degeneracy 
at energy level / is compatible with p{l), the number of integer partitions of /. (e) Schematic representation of the boundary CFT energy 
spectrum of a compactified boson with free boundary conditions [6, 25, 26]. Each orange shaded circle denotes a primary field with scaling 
dimension ri{6NA)^, while the equally spaced white circles on top of each primary field complete the Virasoro towers with a degeneracy count 
of p{l). The schematic spectrum has been plotted for 77 = 1/2. The energy differences denoted A^(0, 0) and A^(0, ±1) in (a) and (c) enter 
the formula (3) ry = A^(0, ±1)/A^(0, 0). 



state Potts model. We then provide a CFT interpration of the 
numerical results and conclude. 

Setup. — We study one-dimensional quantum many body 
systems on chains of length L with open boundary conditions 
(OBC) and periodic boundary conditions (PBC) in the follow- 
ing. The fixed size and position of the block A in the biparti- 
tion is illustrated in Fig. 1 . We calculate the finite size ground 
state wave functions using the density matrix renormalization 
group (DMRG) technique [21] and study the entanglement 
spectrum for various system sizes up to about a thousand sites. 
The results are converged in the number of states kept, where 
this number varies among the models but can be as large as a 
few thousand states. 
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Figure 3. (Color online) Scaling dimension 77 of the Bose-Hubbard 
chain at filling n = N/ L = 1 as a function of the interaction U/J. 
The critical value rjc — 1/4 and location of the transition from the 
superfluid to the Mott insulator {U/ J)c [22, 23] are indicated with 
horizontal and vertical lines. The red dashed curve shows 77 from 
the particle number fluctuations analysis in Ref. [27]. The particle 
number cutoff per site is set to A^max = 4 in all the results. 



Bose Hubbard model. — We investigate the entanglement 
spectrum of Bose-Hubbard chains of length L described by 
the Hamiltonian: 

^BH = -jY.^h\h, + h.c.) + /7/2^n,(n, - 1), (2) 

at fixed unit filling (n^) = 1. Below a critical value ofU/J^ 
3.38 [22] the system is in a ID superfluid phase, described 
by Luttinger liquid theory - i.e. an instance of a c = 1 CFT 
- where some of the scaling dimensions of the theory vary 
continuously with the ratio U/J [23]. 

In Fig. 2(a) we display the ES {^} for a symmetric bipar- 
tition of an OBC wave function at L = 256 and U/J = 2 
using the particle number in the block {SNa '= Na — La) as 
an additional label of the entanglement levels. A prominent 
feature of the ES is that it exhibits a parabolic envelope high- 
lighted by the continuous line through the lowest ES level in 
each SNa sector. In order to uncover the additional structure 
above the parabola we subtract the value of lowest ES in each 
sector, and set the overall scale by assigning the second ES 
level in the SNa = sector a fixed value of one. The result- 
ing spectrum is shown in Fig. 2(b), displaying an intriguing 
equally spaced structure, where furthermore the approximate 
degeneracy at level / seems to be given by the number of 
integer partitions of I. 

Before interpreting these results we show the ES for a chain 
of the same size and Hamiltonian parameters, but with peri- 
odic boundary conditions in Fig. 2(c) [corresponding to setup 
(b) in Fig. 1]. The ES also shows a parabolic envelope, but the 
spectrum is about a factor two denser compared to the OBC 
case. Still, if one applies the shifting and rescaling procedure 
as in the OBC case one finds in Fig. 2(d) the same equally 
spaced spectrum with an identical degeneracy count p{l) at 
level It is quite remarkable that the ES of the PBC setup 
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is isostructural to the OBC case, because in gapped systems 
the structure of the entanglement spectra in the two setups is 
generically different [24] . 

In order to provide an interpretation of the structure of the 
entanglement spectrum it is instructive to recall the form of the 
low energy spectrum of a Bose-Hubbard model in the super- 
fluid phase on a finite system with open boundary conditions. 
For example in Ref. [25] the corresponding energy spectrum 
has been derived using bosonization techniques. The energy 
spectrum has the form sketched in Fig. 2(e), which indeed 
matches the observed structure both in the OBC and PBC en- 
tanglement spectrum. This is one of the key results of this 
paper: the entanglement spectrum has a close correspondence 
to the energy spectrum of a boundary CFT, but with specific 
- here free - boundary conditions, both for OBC and PBC of 
the entire system. 

In the CFT language this spectrum contains a collection of 
primary fields (orange shaded circles) which form the root of 
their respective Virasoro tower of descendent fields [6, 26]. 
Most interestingly the relative position of the primary fields 
expressed in units of the spacing of Lq (denoted here) 
encodes their nontrivial scaling dimensions. Those are re- 
lated to the Luttinger liquid parameter (or in a related way 
the compactification radius of the boson or the decay expo- 
nent of the Green's function). Based on this insight, we can 
now devise a procedure to extract the scaling dimensions di- 
rectly from the ES. If we define the Lq -spacing in the ES by 
A^(0,0) := ii{6NA = 0) - ^o(^A^a = 0) [as indicated in 
Figs. 2(a) and (c)] and the position of the first non-identity pri- 
mary field as A^(0, ±1) := io{SNA = ±1) - io{SNA = 0) 
[as indicated in Figs. 2(a) and (c)], then we can estimate its 
scaling dimension as: 
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In Fig. 3 we plot r] for different system sizes and boundary 
conditions as a function of U/ J. For comparison we also 
plot the estimate for the same quantity obtained from a re- 
cent particle number fluctuations analysis [27]. The agree- 
ment between the different boundary condition choices and 
system sizes and the reference data is very good, particularly 
away from the phase transition. Upon approaching the quan- 
tum phase transition from the superfluid to the Mott insula- 
tor the finite size effects become more pronounced, but still 
theU/J crossing points when r] intersects the r^c = 1/4 line 
(the quantum phase transition is triggered according to Sine- 
Gordon theory when r] crosses r]c = 1/4 [22, 23, 28]) appear 
to converge to the correct transition point as system sizes are 
increased. 

Spin chains. — We now investigate whether the ES based 
extraction of the value of r] is also feasible and accurate for 
quantum spin systems. First we study the 5 = 1/2 XXZ 
model: 
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Figure 4. (Color online) (a) Scaling dimension ?7cc of the 5* = 1 /2 
XXZ model in zero magnetic field (4) obtained from Eq. (3) as a 
function of A = Jz/ Jxy for different system sizes and boundary 
conditions. The dashed line denotes the exact Bethe ansatz result for 
Tjx. (b) Scaling dimension rjx of ih^ S = 1/2 Heisenberg model 
(XXZ model at A = 1) as a function of the magnetization per site. 
The dashed line denotes the numerically exact Bethe ansatz result for 
Tjx from Ref. [31]. (c) The same quantity rjx of the 5* = 1 Heisenberg 
chain as a function of the magnetization per site. The blue circles 
denote the DMRG data of Ref. [34], where the exponent rjx was 
determined based on fits of the decay of the correlation functions. 



This model is integrable by Bethe Ansatz and has a Luttinger 
liquid ground state for — 1 < A = Jz/ Jxy < 1 at = 0. The 
dependence of the scaling dimension rjx on A (governing the 
power-law decay of the transverse spin correlations) is known 
exactly from the Bethe ansatz: r]x{^) = 1 — I/tt arccos(A) 
for = 0. In Fig. 4(a) we show the exact curve together with 
the DMRG ES results for different system sizes and boundary 
conditions. Similar to the Bose-Hubbard case the agreement is 
excellent for small r]x and also improves systematically with 
increasing system size. We tentatively attribute the discrep- 
ancy between the ES estimate for r]x and the exact results at 
larger A to more pronounced finite size effects in the ES due 
to an operator becoming marginal as A 1, similar to the 
Bose-Hubbard chain at the transition to the Mott insulator. Let 
us note that the scheme (3) reproduces the exact result for r]x 
already for finite system size both at A = (r/a^ = 1/2) due 
to the exactly known entanglement spectrum in terms of free 
fermions [29, 30], and at A = 1 (t^^^ = 1), due to the multiplet 
structure of the SU{2) symmetry which is also present in the 
entanglement spectrum. 

As a further application of the method we show the magne- 
tization dependence of r]x for the Heisenberg model (A = 1) 
in a magnetic field in Fig. 4(b). For comparison the numeri- 
cal Bethe ansatz results from Ref. [31] are shown. Finally in 
Fig. 4(c) we present the magnetization dependence of r]x in the 

= 1 Heisenberg chain as a function of the magnetization, 
where an interesting non-monotonous behavior of rjx has been 
reported in past work [32-35]. We also include the results ob- 
tained in Ref. [34] based on fits to the correlation functions 
for comparison, highlighting the good agreement between the 
two complementary approaches. We believe though that our 
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Figure 5. (Color online) Left panel: (a) / (c) shifted and rescaled ES for the transverse field Ising model at the critical point for OBC / PBC. 
The levels which have assigned values are encircled by the dashed red line, (b) catalog of the conformal towers of the Ads c = 1/2 CFT. The 
filled symbols denote the two towers realized in the numerical ES. Right panel: (d) / (f) shifted and rescaled ES for the quantum three- state 
Potts model at the critical point for OBC / PBC. The levels which have assigned values are encircled by the dashed red line, (e) catalog of the 
conformal towers of the A^s c = 4/5 CFT. The filled symbols denote the three distinct towers realized in the numerical ES. 



ES based method is simpler to use than earlier numerical re- 
sults based on correlation function fits. 

In order to ascertain that the correspondence between the 
ES and a boundary CFT extends beyond Luttinger liquids with 
microscopic particle number conservation, we investigate the 
transverse field Ising and the quantum three-state Potts chain 
as further examples. 

The transverse field Ising model. — The transverse field 
Ising chain at its critical point: : 

i i 

realizes the Ising CFT with central charge c = 1/2. DMRG 
results for the entanglement spectrum for both OBC (with free 
EC) and PBC are shown in Fig. 5(a) and (c) respectively. We 
shift the lowest ^ level to zero and subsequently rescale the 
entire spectrum such that the second level is set to 1/2. The 
resulting ES is plotted as a function of 1/ ln(L). The rescal- 
ing is motivated by the structure which becomes apparent in 
the limit 1/ In L ^ 0. The numerical data then plausibly sug- 
gests that the ES arranges on a grid of integer and half integer 
values, with an approximate degeneracy count which can be 
interpreted as the sum of the two Virasoro towers 1/2, see 
Fig. 5(b) for a table of the operator content of the Ising CFT, 
where the realized fields are represented by filled symbols. 
It is remarkable that this identification matches the operator 
content of the partition function Zff of the Ising CFT on an 
annulus with free boundary conditions on both sides [18, 20]. 
We also stress that - as in the compactified boson cases studied 
before - the ES for both OBC and PBC feature the same op- 
erator content. By modifying the boundary conditions in the 
OBC setup (e.g. by applying a magnetic field in the z direc- 
tion at both ends) it is possible to change the operator content 
and reveal a single tower with scaling dimension 1/16 [36], in 



agreement with [18]. 

Finally we note that one can also calculate the ES using the 
mapping to free fermions [30, 37], that a relation between Vi- 
rasoro characters and the spectrum of reduced density matri- 
ces based on the corner transfer matrix off criticality has been 
discussed previously [38, 39] and in Ref. [40] the spectrum 
of a critical corner transfer matrix was studied in a different 
setup, but a counting similar to ours was observed. 

The quantum three- state Potts model. — As a last example 
we study the quantum three-state Potts model on a chain at 
its critical point [41]. This lattice model realizes a CFT with 
central charge c = 4/5. We plot the ES data obtained by 
DMRG using OBC (with free BC) and PBC based on a sim- 
ilar rescaling procedure, but this time we choose to scale the 
fourth level to a fixed value of two. The resulting ES spectra 
are displayed in Fig. 5(d) and (f). This rescaling is justified a 
posteriori by the emergent structure. Comparing the catalog 
of primary fields and their descendents of the minimal 
model CFT [6, 7] [shown in Fig. 5(e)] with the ES data it is 
plausible that the ES has the operator content: 0© (2 x 2/3) ©3. 
This operator content is identical to the one of the CFT annu- 
lus partition function Zff of the Potts model with free bound- 
ary conditions [18, 19]. By modifying the boundary condi- 
tions in the OBC setup (e.g. by applying a local field prefer- 
ring or suppressing one of the three states at both ends) it is 
possible to change the operator content and reveal the towers 
1/8 © 13/8 or 1/40 © 21/40 [36], in agreement with [18, 19]. 

CFT interpretation. — We have provided substantial numer- 
ical evidence that the entanglement spectrum of quantum lat- 
tice models at a conformal critical point can be interpreted 
as the spectrum of the Lq operator in a boundary CFT. The 
idea that the ES of a region is related to Lq appears repeat- 
edly in the CFT entanglement literature [8, 39, 42-44], but 
the specific boundary conditions and the resulting operator 



5 



content have not been clearly stated. In the models studied 
in the present work we find the conformal boundary condi- 
tions to be free. The identical operator content of the OBC 
and PBC setup provides further evidence that for these sys- 
tems the boundary condition at the entanglement cut is a free 
boundary condition. Free boundaries at the entanglement cut 
also appear naturally in the replica approach by Calabrese and 
Cardy [11, 45]. An analytical derivation for the numerical ob- 
servations reported here is needed, and might also address the 
origin and possible additional useful information encoded in 
the finite- size corrections to the ES (such as marginal opera- 
tors). It will also be important to understand whether differ- 
ent microscopic implementations of the same CFT can lead to 
different boundary conditions in the boundary CFT and thus 
modify the operator content of the ES. 

In this paper we have shown that the real space entangle- 
ment spectrum gives access to the operator content (or equiv- 
alently the energy spectrum) of certain boundary CFTs. In 
the particular case of those Luttinger liquids considered in 
this work it is possible to extract the Luttinger liquid param- 
eter solely based on the entanglement spectrum. This tech- 
nique has the potential to become a standard diagnostic tool in 
DMRG or related matrix product state based techniques (such 
as iTEBD [46], iDMRG [47] or cMPS [48], possibly with dif- 
ferent CFT boundary conditions and operator content [36]) to 
characterize one-dimensional conformally invariant systems 
or multicomponent generalizations thereof (such as fermionic 
Hubbard models with gapless spin and charge degrees of free- 
dom). 
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